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Abstract:  Molecular  orbital  calculations  of  the  chemisorption 
and  reactivity  of  ethylene  sulfide  and  trimethylene  sulfide  on 
Mo (110)  are  presented,  and  compared  to  similar  binding  in  model 
discrete  Mo  complexes.  Our  calculations  suggest  preferred 
bonding  of  the  cyclic  sulfides  on  2-  or  3-fold  sites  of  the 
surface,  by  the  expected  S-lone  pair  donor  mechanism.  The 
concerted  elimination  of  ethylene  or  cyclopropane  is  much  easier 
on  the  surface  than  it  is  in  model  organometallic  molecules.  The 
activation  barrier  for  ethylene  sulfide  decomposition  calculated 
is  substantially  smaller  than  for  trimethylene  sulfide.  Various 
nonconcerted  mechanisms  and  the  role  of  coadsorbed  sulfur  and 
hydrogen  were  also  probed. 
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Hydrodesulfurization  is  a  widely  used  process,  by  which  fuel 
feedstocks  react  with  a  molybdenum  sulfide  catalyst  to  afford 
hydrocarbons  with  lower  sulfur  content.  Several  kinds  of  sulfur 
derivatives,  both  aromatic  and  aliphatic,  are  found  in  the  raw 
materials  and  considerable  effort  has  been  dedicated  towards 
understanding  how  the  reactions  occur.1  Often,  simpler  models 
than  the  actual  catalyst  are  used  for  that,  namely  single  crystal 
transition  metal  surfaces.2 

The  molybdenum  (110)  face  has  been  chosen,  due  to  its 
stability  towards  reconstruction,  by  Friend  and  coworkers  to 
study  systematically  the  reactions  of  cyclic  sulfides3  and  linear 
thiols.*  With  the  help  of  several  complementary  experimental 
techniques  such  as  temperature  programmed  desorption.  X-ray 
photoelectron  and  high  resolution  electron  energy  loss 
spectroscopy,  they  have  been  able  to  establish  three  distinct 
reactivity  patterns.  These  three  pathways  are  indicated  in  1, 


A 
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depicting  how  trimethylene  sulfide  reacts  after  adsorption  onto 
the  Mo (110)  face. 

Pathway  A  leads  to  cyclopropane  formation  by  intramolecular 
elimination.  An  intermediate  adsorbed  thiolate  species  is  formed 
in  B.  The  same  intermediate  is  detected  directly  when  a 
different  precursor,  propanethiol ,  reacts  with  the  surface. 

The  adsorbed  thiolate  species  can  react  further  breaking  the 
C-S  bond,  either  by  forming  another  C-H  bond,  or  losing  another  H 
to  the  surface,  yielding  respectively  alkanes  or  alkenes. 

Complete  decomposition  occurs  in  pathway  C,  resulting  in  surface 
sulfur,  carbon,  and  hydrocarbon  fragments. 

This  kind  of  reactivity  appears  to  be  quite  general  for 
adsorbed  thiolates:  for  instance,  following  the  same  trends, 
phenyl thiolate  adsorption  results  in  gaseous  benzene  and  surface 
benzyne.  However,  cyclic  sulfides,  depending  on  ring  size,  react 
differently  on  Mo(110) .  Pathway  B  is  not  observed  for  ethylene 
sulfide,  whereas  intramolecular  elimination  (A)  does  not 
occur  for  ring  sizes  greater  than  four. 

We  would  like  to  understand  how  the  different  sulfur 
containing  species  interact  with  the  molybdenum  surface  and  what 
electronic  factors  govern  their  reactivity.  Let  us  note  that 
these  reactions  occur  in  the  presence  of  hydrogen,  which  is  both 
an  impurity  impossible  to  eliminate,  and  a  product  from  S-H  bond 
cleavage  in  the  case  of  thiol  reactions.  Also,  as  reactions 
proceed,  more  and  more  surface  carbon  and  sulfur  are  formed, 
changing  the  reactive  capabilities  of  the  surface. 3f 
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Gaining  some  knowledge  about  the  reactivity  of  sulfur- 
containing  species  on  a  molybdenum  (110)  surface  may  help  us  to 
understand  how  real  catalytic  hydrodesulfurization  takes  place 
under  more  complicated  experimental  conditions.  Many  of  the 
studies  done  concern  thiophene,  which  is  considered  one  of  the 
most  difficult  molecules  to  desulfurize.2,3  Molecular  thiophene 
complexes  have  been  synthesized  and  used  as  models  for 
hydrodesulfurization  behavior,  with  variable  success.6 

Some  discrete  molecular  complexes  are  available,  and  we 
shall  use  them  as  models  for  coordination  sites  and  geometries. 
Also,  understanding  how  a  cyclic  sulfide  or  thiolate  interacts 
with  a  mono-or  binuclear  molybdenum  complex  will  help  us  to 
interpret  its  interaction  with  a  two-dimensional  surface.  These 
two  approaches,  molecular  and  extended  surface,  are  conceptually 
different,  but  should  be  complementary  in  their  results. 

The  calculations  done  are  of  the  extended  Hiickel  type,7  and 
the  tight-binding  method8  was  used  to  study  the  two-dimensional 
infinite  surface.  Details  are  given  in  the  text  and  in  the 
Appendix. 
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The  adsorbates 

In  the  following  study,  we  shall  consider  two  adsorbates, 
SC2H4,  ethylene  sulfide,  and  SC3H6,  trimethylene  sulfide,  as 
models  for  the  different  reactivity  patterns  of  cyclic  sulfides 
on  Mo ( 110) . 

So  far  nothing  is  known  about  the  geometry  of  these 
heterocycles  after  adsorption  on  the  surface.  One  discrete 
molecular  complex  of  ethylene  sulfide  has  been  characterized  by 
Rauchfuss  and  coworkers.9  It  is  a  CpRu(PPh3)2(SC2H4)  +  molecule, 
and  the  geometry  of  the  ethylene  sulfide  ring  in  it  is  similar  to 
that  of  the  free  molecule.  In  our  calculations  we  take  the 
experimentally  observed  free  molecule  structures,  from  microwave 
spectroscopy  for  ethylene  sulfide10*  and  from  gas  electron 
diffraction  for  trimethylene  sulfide. 10b 

Ethylene  sulfide  is  a  planar  molecule  with  S-C  bonds  1.815A 
long  and  a  C-C  bond  of  1.484A.  Its  frontier  molecular  orbitals 
are  shown  in  2.  ’The  LUMO,  2blf  both  S-C  and  C-C  antibonding,  is 
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quite  high  in  energy.  The  HOMO,  lb2,  and  SHOMO,  2alf  are 
respectively  the  n  and  a  sulfur  lone  pairs,  while  the  next 
orbital,  lblf  is  strongly  S-C  bonding  but  C-C  antibonding. 
Depopulation  of  this  orbital  may  be  important  in  weakening  the  S- 
C  bond,  and  simultaneously  strengthening  the  C-C  bond.  This 
would  put  the  system  on  the  way  to  produce  gaseous  ethylene  and 
surface  sulfur. 

Trimethylene  sulfide  was  experimentally  found  to  be  non- 
planar.  The  folding  angle  between  the  C-C-C  and  C-S-C  planes  is 
26*,  the  S-C  and  C-C  bond  lengths  respectively  1.85  and  1.549A. 
The  two  sulfur-bonded  carbon  atoms  are  2.294A  apart.  The  planar 
four  membered  ring  is  about  0.6eV  higher  in  energy,  according  to 
our  calculations. 

The  frontier  molecular  orbitals  for  folded  trimethylene 
sulfide  are  given  in  3.  The  LUMO  lies  at  a  higher  energy  than 
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for  ethylene  sulfide  and  the  two  highest  occupied  M.O.'s  are 
again  the  ir  and  a  sulfur  lone  pairs. 

These  molecules  will  bond  to  a  metal  atom  essentially 
through  the  sulfur  lone  pairs,  but  lower  molecular  orbitals  may 
be  involved  in  some  situations.  Let  us  consider  ethylene  sulfide 
and  two  limiting  coordination  geometries,  shown  in  4. 


The  total  energy  changes  while  going  from  the  perpendicular 
coordination  (4,  left)  to  the  parallel  one  (4,  right).  This  is 
shown  in  5,  for  MLn=Mo(CO)s,  MoH*4*’  and  CH3+. 


180  160  140  120  100 
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In  general  sulfonium  salts  are  strongly  pyramidal,  and  our 
calculations  reproduce  this  fact.  Steric  constraints  play  a 
role,  as  the  ethylene  sulfide  hydrogens  approach  other  atoms. 

This  is  why  the  small  CH3  group  gives  the  largest  bending  angle. 
In  the  one  known  organometallic  complex  of  ethylene  sulfide,  a  is 
112*  .9 

The  way  SC2H4  binds  to  a  square  planar  MoH*4'  (d6,Mo(0)  ,  in 
anticipation  of  the  Mo  surface)  is  analyzed  in  the  interaction 
diagram  of  Figure  1  for  the  perpendicular,  higher  symmetry, 
geometry.  On  the  right  side  are  the  orbitals  of  SC2H4,  as  shown 
before  in  2.  The  orbitals  on  the  left  are  those  of  a  square 
planar  complex,  with  the  typical  four  below  one  splitting  of  the 
d  levels.11  The  highest  one  is  not  seen  in  this  energy  window. 

Figure  l  here 

Large  overlaps  are  present  between  the  sulfur  a  lone  pair 
(2ai  and  lai)  and  empty  3ax  (z2)  and  4a!  (z)  of  the  metallic 
fragment.  Back  donation  from  2bj  (xz)  to  the  ethylene  sulfide 
LUMO  is  negligible.  The  sulfur  ligand  acts  as  a  a  donor.  The 
HOMO  is  involved  in  a  weak  four-electron  destabilizing 
interaction. 

When  the  three-membered  ring  bends,  the  symmetry  is  lowered 
and  the  HOMO  can  interact  with  3a:  and  4ai  of  MoH/'  There  is  a 
better  energy  matching  between  these  orbitals  and  a  stronger 
interaction  results.  The  previous  destabilizing  interaction  is 
relieved.  A  better  bond  is  formed,  as  shown  by  the  minimum  in 
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Figure  1.  Interaction  diagram  for  the  model  complex 
MoH*  (SC2H*) 
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the  total  energy  curve  (5)  for  large  bending  angles. 

Let  us  now  consider  ethylene  sulfide  acting  as  a  bridging 
ligand.  There  are  two  extreme  geometries  —  parallel,  6,  and 
perpendicular,  7.  Furthermore  one  could  imagine  bending  starting 
from  either  conformation. 
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Although  in  the  parallel  form  some  stabilization  is  achieved 
by  bending,  the  initial  perpendicular  geometry  is  much  more 
stable.  In  comparing  the  perpendicular  and  parallel  geometries, 
the  a  interaction  between  the  2at  lone  pair  and  symmetric 
combinations  of  Mo  z2  and  xz  orbitals  will  be  similar.  However, 
for  the  perpendicular  geometry  a  strong  interaction  develops,  8, 
between  the  antisymmetric  combinations  of  z2  and  xz  and  the  HOMO. 
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An  analogous  interaction  for  the  parallel  geometry  would  involve 
lbj,  which  is  a  poorer  donor  and  has  a  smaller  localization  on 
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the  sulfur.  It  would  be  weaker  due  both  to  energy  and  overlap 
reasons.  The  resulting  species  would  then  be  less  stable. 

An  example  of  a  bridging  cyclic  sulfide  (tetra- 
hydrothiophene)  adopting  a  perpendicular  geometry  is  seen  in  the 
complex  [Nb2Cl3  (ju-tetrahydrothiophene)  3]2'1  9. 12 


2- 
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Our  calculations,  not  reported  here  in  detail,  show  that  the 
bonding  preferences  of  trimethylene  sulfide,  SC3H6,  are  similar 
to  those  of  ethylene  sulfide. 

The  molecular  complexes  used  as  models  are  related  to  the 
molybdenum  environment  in  the  surface.  The  symmetry  is  higher  in 
these  discrete  molecules,  making  the  analysis  simpler,  but  the 
results  are  qualitatively  the  same  compared  to  more  realistic 


models. 
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The  reactivity  of  the  adsorbates 

Let  us  imagine  a  concerted  mechanism  which  leads  from  an 
ethylene  sulfide  complex  to  a  sulfide,  with  formation  of  free 
ethylene,  10,  in  a  process  similar  to  that  proposed  on  the 
Mo (110)  surface. 
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Similarly,  we  shall  have  11  for  trimethylene  sulfide, 
leading  to  cyclopropane. 


II 

The  related  reaction  between  ethylene  oxide  and  the  tungsten 
complex  WCl2(PMePh2K  has  been  studied  theoretically.13  However, 
this  reaction  differs  from  the  one  we  want  to  study  in  that, 
after  breaking  the  C-0  bonds,  both  the  oxygen  and  the  ethylene 
remain  bound  to  the  metal  center.  In  the  surface  reaction,  no 
ethylene  is  evolved  into  the  gas  phase  at  the  adsorption 
temperature  and  none  is  bound  directly  to  the  surface,  suggesting 
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that  the  gaseous  product  does  not  proceed  through  a  surface-bound 
ethylene  intermediate.  The  reaction  dynamics  are  suggestive  of  a 
concerted  transition  state. 

Let  us  then  assume  a  concerted  reaction:  the  two  S-C  bonds 
stretch  simultaneously,  while  the  hydrogens  relax  towards  a 
planar  ethylene  geometry  by  the  end  of  the  reaction,  as  indicated 
in  12.  The  C-C  bond  length  was  kept  constant. 


ML„ 
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contracting  it  from  1.484A  to  1.34A  (in  free  ethylene)  will  lead 
to  further  stabilization.  The  change  of  total  energy  and  the 
Walsh  diagram  along  this  hypothetical  reaction  path  are  shown  in 
Figure  2,  as  a  function  of  the  S-C  bond  stretching. 

Figure  2  here 

C2v  symmetry  is  retained  throughout  the  reaction.  There  are 
several  avoided  crossings.  One  originates  from  a  high  energy  s-c 
antibonding  orbital  that  cannot  correlate  with  the  n  orbital  of 
the  forming  ethylene,  which  is  much  more  stable.  Initially, 
there  are  six  electrons  in  the  three  lowest  d  orbitals.  By  the 
end  of  the  reaction,  n  will  be  occupied  and  only  four  electrons 


H  — Mo'— H 
YT 


H— Mo— H  +• 


Figure  2 


Walsh  diagram  and  total  energy  change  for  the  reaction 
MoH*  (SC2HJ  *'  -*  M0H4S*"  +  C2Ht,  as  a  function  of  the 
stretching  of  the  S-C  bond. 
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remain  in  the  metal  d  orbitals.  There  is  a  formal  oxidation 
state  change  from  Mo(0)  to  Mo(II),  as  the  sulfide  complex  is 
formed.  The  reaction  is  forbidden  for  this  electron  count,  with 
a  computed  activation  energy  close  to  lev. 

Let  us  compare  these  results  to  those  obtained  when  the 
reacting  ethylene  sulfide  is  bridging  two  molybdenum  atoms,  in 
the  model  compound  Mo2H8(SC2HJ 6"  (perpendicular  geometry),  7. 

The  metal  fragment  is  derived  from  two  square  planar  units 
and  eight  d  orbitals  are  now  present,  the  other  two  having  too 
high  energies.  The  symmetry  is  still  C2v  and  the  assumed 
reaction  pathway  is  exactly  the  same  as  before.  The  Walsh 
diagram  and  the  change  in  total  energy  are  shown  in  Figure  3 . 
Generally,  the  energy  varies  in  the  same  way  as  was  observed  in 
the  previous  model  and  analogous  level  crossings  are  observed 
here  between  orbitals  of  ax  symmetry. 

Figure  3  here 

The  same  essentially  holds  if  we  take  other  models  that  more 
closely  resemble  the  type  of  molybdenum  coordination  on  the  (110) 
face,  with  six  nearest  neighbors,  four  in  the  surface  plane  (but 
not  arranged  in  a  square)  and  two  below. 

What  happens  if  we  now  look  at  the  reaction  represented  in 
11,  that  is,  formation  of  cyclopropane  and  a  sulfide  complex  from 
the  trimethylene  sulfide  complex?  The  model  is  MoH*(SC3H6) and 
the  trimethylene  sulfide  ligand  is  bound  as  shown  in  the  right 
side  of  4,  with  a  M-SCC  plane  angle  of  130*. 


Energy  (eV) 


Figure 


.  Walsh  diagram  and  total  energy  change  for  the  reaction 
Mo2H8(SC2HJ  6‘  -*  Mo2H6S6‘ +  C2H*,  as  a  function  of  the 
stretching  of  the  S-C  bond. 
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We  must  start  by  defining  a  reaction  path.  Here  there  are 
more  severe  geometrical  changes  to  contend  with.  While  for 
ethylene  sulfide  a  single  carbon-carbon  bond  became  a  double 
bond,  now  a  new  carbon-carbon  bond  has  to  be  formed  starting  from 
atoms  that  are  relatively  far  apart.  In  13  are  shown  the 
geometrical  parameters  to  be  varied:  the  S-C  bond  is  stretched; 
the  two  carbon  atoms  approach  each  other  to  form  a  bond;  the 
hydrogens  bound  to  them  reorient. 
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In  the  end  the  three  C-C  bonds  in  the  cyclopropane  product 
have  the  same  length  (but,  as  for  ethylene,  we  did  not  allow  for 
relaxation  to  the  true  C-C  bond  length) .  The  symmetry  is  now 
lower  than  for  the  formation  of  ethylene.  Due  to  the  folding  of 
the  SC3H6  ring,  only  one  mirror  plane  is  left  to  the  molecule, 
instead  of  two. 

The  Walsh  diagram  is  very  similar  to  the  one  in  Figure  2 . 
Instead  of  a  n  orbital,  a  stabilized  level  corresponding  to  the 
newly  formed  C-C  bond  is  observed  at  low  energies  towards  the  end 
of  the  reaction.  It  derives  from  a  high  energy  antibonding 
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orbital,  but  avoided  crossings  occur  with  metal  d  orbitals  of  the 
same  symmetry.  The  shape  of  the  total  energy  curve,  14,  is 
similar  to  the  ones  seen  before.  The  reaction  is  still 
forbidden,  with  a  relatively  small  activation  energy. 


Let  us  now  look  at  the  evolution  of  the  S-C,  Mo-S  and  C-C 
overlap  populations  along  the  reaction  pathway,  as  they  give  us 
an  indication  of  how  bonds  are  being  formed  or  broken. 

They  are  represented  in  Figure  4  for  the  two  reacting 
species:  ethylene  sulfide  in  Mo2H6  (SC2H6) 6'  and  trimethylene 
sulfide  in  MoH* (SC3H6) 4*  There  are  no  great  changes  in  the  Mo-S 

Figure  4  here 

bond  while  the  S-C  bond  is  broken  and  the  C-C  is  strengthened 
(for  SC2HJor  formed  (SC3H6) .  The  changes  are  relatively  smooth, 
some  discontinuities  arising  from  the  above  mentioned  avoided 
crossings.  The  proposed  pathway  seems  reasonable. 

We  are  now  ready  to  proceed  to  the  study  of  the  surface 


reactions. 


Figure  4.  Evolution  of  the  S-C,  C-C,  and  Mo-S  overlap  populations 
along  the  pathways  described  in  the  text  for  the 
reactions  Mo2H6(SC2HJ 6' -  Mo2H6S6‘ +  C2H*  (left)  and 
MoH*(SC3H„)*‘ -  MoHtS6"  +  C3Hs  (right). 


16 


The  Mo  ( llCH  Surface 

Molybdenum  has  a  body  centered  cubic  structure  with 
a=3.15A.u  The  110  face  is  centered  rectangular,  but  we  use  a 
P (2x2 )  unit  cell,  IS. 


3  15  l 
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The  atoms  in  the  second  layer  occupy  positions  in  the  center 
of  the  rhombus,  16,  and  the  stacking  is  ABAB. . . ,  that  is  the 
third  layer  is  exactly  below  the  first  one. 


•  X  •  X  • 

X  •  X  •  X 

•  X  •  X  • 

X  •  X  •  X 

•  X  •  X  • 
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We  model  the  surface  by  a  slab  of  three  layers  of  molybdenum 
atoms.  This  number  has  been  shown  to  be  a  reasonable  compromise 
between  a  convenient  size  of  the  calculation  and  an  adequate 
representation  of  the  surface.15 

The  choice  of  the  unit  cell  was  dictated  by  the  necessity  of 
keeping  a  relatively  low  coverage  on  the  surface,  otherwise  the 
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adsorbate  molecules  may  come  into  close  contact  with  one  another, 
leading  to  high  repulsive  energies.  With  the  chosen  unit  cell, 
the  coverage  will  be  1/4  (one  adsorbate  molecule/ four  surface 
molybdenum  atoms)  and  no  short  distances  between  adjacent 
adsorbate  molecules  are  observed. 

We  took  the  Mo-S  distance  to  be  2.44A,  which  is  a  very 
typical  one,  found  in  many  molybdenum  complexes.3*1'16  This 
constant  distance  implies  that  the  height  of  the  coordinating 
sulfur  atom  in  the  adsorbate  above  the  surface  will  vary 
according  to  the  multiplicity  of  the  site.  In  what  follows,  we 
are  going  to  consider  three  different  adsorption  sites:  on  top, 
17a;  bridging  between  adjacent  molybdenum  atoms  (nearest 
neighbors) ,  17b;  3-fold,  with  the  sulfur  atom  equidistant  from 
Mol,  Mo2,  and  Mo4,  17c.  Having  picked  the  adsorption  sites,  we 

a 


b 


c 

y  |  *■»  Adsorbate 
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have  to  decide  on  the  adsorption  geometry.  Both  ethylene  sulfide 
and  trimethylene  sulfide  contain  a  planar  SCC  unit.  From  what  we 
saw  about  their  frontier  orbitals,  2  and  3,  we  may  expect 
different  interactions  with  one  (or  more)  molybdenum  atoms 
depending  on  the  orientation  of  that  ring.  We  shall  consider, 
both  for  the  top  and  3-fold  sites,  limiting  adsorbate  geometries 
with  the  ring  in  the  xz  or  yz  plane;  the  corresponding  situation 
for  the  2-fold  site  will  have  the  ring  in  the  same  plane  as  the 
Mo-Mo  bond  or  perpendicular  to  it.  These  possible  situations  are 
sketched  in  18. 


’L 
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Allowing  the  ring  to  bend  introduces  other  interesting 
possibilities,  which  will  be  discussed  later. 
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Adsorption  of  ethvlene  sulfide 

A  first  step  in  trying  to  understand  how  ethylene  sulfide 
reacts  on  the  molybdenum  (110)  surface  involves  finding  its 
preferred  adsorption  site.  In  order  to  do  that,  we  look  at 
different  parameters,  such  as,  for  instance,  the  binding  energy 
of  the  adsorbate  to  the  surface  and  the  overlap  population 
between  the  sulfur  atom  of  the  adsorbate  and  the  molybdenum  to 
which  it  binds.  These  and  other  important  parameters  are 
collected  in  Table  1  for  several  possible  situations.  A  "  bent" 
adsorbate 

Table  1  here 

means  that  the  two  carbons  and  their  hydrogens  move  towards  the 
surface,  as  measured  by  an  angle,  a,  from  the  vertical,  19. 
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Let  us  start  analyzing  the  results  by  looking  at  the  binding 
energies.  These  are  all  positive,  meaning  that  a  stabilizing 
interaction  has  taken  place  between  the  adsorbate  and  the 
surface.  Notice  that  the  smaller  values  occur  for  on-top 
adsorption,  where  the  adsorbate  prefers  to  be  bent.  The  opposite 
happens  in  the  3-fold  site,  where  bending  is  decidedly 


Binding  energies  and  overlap  populations  for  thiirane  on  Mo (110) 


(binding  energy)  =  E(Mo  slab)  +  E(thiirane  layer)  -  E(Mo  +  thiirane) 


20 


unfavorable.  This  is  the  same  trend  that  we  observed  for  the 
bonding  of  ethylene  sulfide  to  mono-and  binuclear  molybdenum 
complexes.  The  binding  energies  for  the  2-fold  and  3-fold  sites 
are  large  and  we  cannot  trust  the  difference  between  them,  since 
they  will  depend  on  the  exact  Mo-S  bond  lengths,  etc.  We  can 
conclude  that  all  these  sites  are  equally  probable,  more  so  than 
the  top  site.  The  bent  3-fold  site  is  also  probably  excluded. 

Strong  molybdenum  sulfur  bonds  are  formed  in  all  cases. 
Although  the  largest  Mo-S  overlap  population  is  found  for  the  top 
site,  only  one  bond  per  unit  cell  is  formed.  For  the  2-fold 
site,  two  slightly  weaker  bonds  are  formed  and  three  weaker  ones 
are  found  for  the  different  possibilities  of  3-fold  adsorption. 
Here,  two  different  numbers  are  shown,  as  the  bonds  from  S  to  Mo2 
and  Mo4  are  different  from  the  bond  between  Mo  and  S,  even  if  the 
Mo-S  distance  is  the  same.  While  it  is  difficult  to  draw 
conclusions  from  these  overlap  populations,  it  seems  that  a 
stronger  bond  to  the  surface  is  formed  in  the  2-fold  and  3-fold 
adsorption  sites. 

The  remaining  overlap  populations  in  Table  1  give  us  some 
information  about  what  happened  to  the  ethylene  sulfide  molecule 
after  adsorption.  Changes  in  bond  strengths  due  to  the 
interaction  with  the  surface  may  be  indicative  of  the  new 
reactivity  the  molecule  exhibits.  If  we  consider  the  free  three- 
membered  ring,  the  S-C  overlap  population  is  0.560  and  the  C-C 
overlap  population  0.737.  For  both  top  and  2-fold  sites,  the  S-C 
overlap  population  remains  unchanged,  while  the  C-C  weakens. 
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This  site  is  probably  not  a  precursor  to  gaseous  ethylene 
formation.  On  the  other  hand,  a  simultaneous  weakening  of  the  S- 
C  bond  and  strengthening  of  the  C-C  is  observed  when  ethylene 
sulfide  binds  to  the  surface  in  the  3-fold  site,  suggesting  an 
activation  towards  formation  of  ethylene,  as  is  experimentally 
observed. 

In  order  to  understand  these  results,  let  us  analyze  the 
bonding  of  one  of  these  sites  in  more  detail.  We  can  choose,  for 
instance,  the  situation  where  the  three-membered  ring  adsorbs  in 
the  3 -fold  site,  perpendicular  to  the  surface,  with  its  ring 
lying  parallel  to  yz. 

We  show  in  the  center  of  Figure  5  the  total  density  of 
states  (DOS)  of  a  Mo (110)  slab  covered  by  an  adsorbate  layer. 

The  dark  area  gives  a  projection  of  all  the  adsorbate  states. 
These  states  were  obtained  by  interaction  of  the  Mo (110)  slab 
states,  shown  on  the  left  side  of  Figure  5,  with  those  of  the 

Figure  5  here 

adsorbate  layer,  given  on  the  right  side.  The  molybdenum  DOS  is 
characteristic  of  its  metallic  properties,  with  no  gaps.  The  s 
band  penetrates  the  d  band  which  lies  between  -8  and  -13eV. 

Sharp  peaks  are  observed  for  the  adsorbate  layer  DOS,  as  each 
molecule  is  quite  far  from  its  neighbors,  thus  resulting  in  very 
weak  interactions.  These  sharp  peaks  can  develop  into  broader 
bands  by  interaction  with  the  Mo  surface.  Each  peak  is  labelled 
according  to  the  molecular  orbital  it  is  derived  from  in  2. 


(A*)  X8j«u3 
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Comparing  the  three  parts  of  Figure  5,  we  can  see  more 
easily  changes  in  the  adsorbate  DOS  than  in  the  molybdenum  layer 
DOS.  These  indicate  that  some  interaction  is  taking  place. 
However,  to  understand  which  orbitals  are  involved  in  the  bonding 
we  must  consider  the  projections  of  each  of  them  individually. 

This  is  shown  in  Figure  6  for  several  important  ethylene 
sulfide  orbitals.  In  (a)  we  see  the  projection  of  the  three- 

Figure  6  here 

membered  ring  HOMO,  lb2,  solid  line,  with  the  respective 
integration  (dotted  line) .  There  is  still  one  main  peak,  but  the 
integration  shows  this  fragment  orbital  to  be  spread  over  a  large 
energy  range;  -15%  of  the  states  are  localized  in  the  Mo  d 
region.  In  addition,  the  main  peak  is  pushed  down  from  ~l2.5eV 
to  ~l3.5eV,  due  to  interaction  with  the  surface.  lb2  interacts 
especially  effectively  in  the  three-fold  site.  Something  similar 
is  happening  with  the  other  two  orbitals  shown  in  Figure  6  (b) 
and  (c) ,  which  are  respectively  the  other  sulfur  lone  pair  2ax 
and  the  molecular  orbital  lax.  After  adsorption,  some  of  the 
ethylene  sulfide  states  are  partially  vacated,  their  electrons 
transferred  to  molybdenum  states.  This  transfer  of  electrons 
from  the  sulfur  ligand  to  the  metal  can  also  be  seen  in  the 
occupation  of  the  levels  mentioned,  shown  in  Table  2,  which 
includes  the  other  adsorption  sites  as  well.  We  can  see  that 
the  lbx  level  is  only  significantly  depopulated  in  the  3-fold 
sites.  2ax  and  lb2  are  depopulated  for  all  of  the  sites, 


(A*)  «J*ug 


three  ethylene  sulfide  orbitals  after  interaction 
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although  lb2  is  more  so  for  the  higher  coordination  sites. 


Table  2  here 

The  bonding  characteristics  of  the  adsorbate  levels  are 
shown  in  the  Crystal  Orbital  Overlap  Population  (COOP)  curve  of 
Figure  7.  While  the  HOMO  (lb2)  is  essentially  nonbonding  and  2al 
is  C-C  bonding,  lbx  is  simultaneously  S-C  bonding  and  C-C 
antibonding.  If  it  is  depopulated,  the  S-C  bond  will  be  weaker 
and  the  C-C  bond  stronger.  That  is  a  step  towards  breaking  the 
S-C  bonds  and  producing  ethylene. 

Figure  7  here 

The  ethylene  sulfide  LUMO  (2^)  interacts  strongly  with 
surface  levels,  as  can  be  seen  in  Figure  5  (the  peak  of  5c  is 
completely  broken  up) .  However  this  level  is  only  weakly 
populated,  as  the  states  it  gives  rise  to  are  essentially  above 
the  Fermi  level.  This  is  important  because  the  major  cause  of 
weakening  of  C-S  and  strengthening  of  C-C  bonds  upon  adsorption 
is  not  population  on  this  LUMO,  but  depopulation  of  lb^ 


Table  2 


Occupation  of  thiirane  orbitals  after  adsorption  on  Mo (110) 


1 

Ms3 

N 

S 

-b 

y 

■b 

<1 

-b 

2b! 

(LUMO) 

0.040 

0.034 

0.077 

0.059 

0.052 

0.063 

lb2 

(HOMO) 

1.885 

1.815 

1.594 

1.538 

1.552 

1.616 

2al 

1.839 

1.886 

1.829 

1.830 

1.829 

1.817 

lbi 

1.968 

1.969 

1.944 

1.806 

1.744 

1.785 
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Figure  7.  Crystal  Orbital  Overlap  Population  (COOP)  curve 
depicting  the  bonding  properties  of  the  ethylene 
sulfide  levels  relative  to  the  S-C  and  C-C  bonds. 
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Adsorption  of  trimethvlene  sulfide 

Trimethylene  sulfide  is  not  expected  to  behave  very 
differently  from  ethylene  sulfide  as  both  its  geometry  (around 
the  coordinating  sulfur  atom)  and  its  frontier  orbitals  (3)  are 
similar.  Some  results  concerning  adsorption  in  2-  and  3-fold 
adsorption  sites  are  collected  in  Table  3. 

Table  3  here 

The  binding  energies  are  positive  both  for  the  2-fold  and  3- 
fold  sites,  but  they  become  much  smaller  when  the  adsorbate  is 
bent  50*  down  from  perpendicular.  Strong  Mo-S  bonds  are  formed 
(more  nonequivalent  Mo-S  bonds  are  now  formed,  due  to  the 
nonplanarity  of  the  adsorbate) ,  in  a  comparable  way  to  ethylene 
sulfide.  The  S-C  bond  is  weakened  by  adsorption  on  the  3-fold 
site,  but  not  on  the  2-fold,  as  can  be  seen  from  the  respective 
overlap  populations.  The  C-C  bonds  are  kept  mostly  unchanged, 
except  for  the  3-fold  adsorption  with  the  C-S-C  plane  parallel  to 
yz.  If  cyclopropane  is  to  be  formed  following  adsorption  the  two 
a-carbons  have  to  approach  each  other  to  form  a  bond.  The 
overlap  population  between  them,  denoted  as  C*"C,  is  always 
negative,  although  small.  Adsorption  does  not  seem  to  activate 
this  molecule,  as  opposed  to  what  we  saw  happening  for  ethylene 
sulfide. 

In  Figure  8  are  shown  both  the  adsorbate  layer  (right)  and 
molybdenum  slab  DOS  (left)  before  adsorption,  and  the  total  DOS 
for  the  composite  structure  (center) .  Significant  changes  take 


Table  3 


Binding  energies  and  overlap  populations  for  thietane  on  Mo(lio) 


Adsorption  site, 

geometry  18d,  perp.  18f,  perp.  18e,  perp.  18e,  bent  505 


B.  E.  (eV)  2.66  2.30  2.75  1.75 


Overlap  population 

S-C  0.628 

(free  SC3H6=0.627) 

C-C  0.705 

(free  SC3H6=0.708) 

C- • •  C  -0.099 

(free  SC3H6=-0 . 099) 

Mol-S  0.518 

Mo2-S  0.518 

Mo3-S 


0.603 

0.608 

0.597 

0.693 

0.707 

0.708 

0.089 

-0.087 

-0 . 086 

0.469 

0.322 

0.320 

0.366 

0.457 

0.353 

0.366 

0.463 

0.450 
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place,  indicative  of  the  interaction  between  the  two  layers.  A 


Figure  8  here 

more  detailed  look  at  the  adsorbate  orbitals  reveals  that 
interaction  proceeds  mainly  through  the  two  highest  occupied 
molecular  orbitals,  as  they  suffer  the  largest  depopulation  and 
the  states  spread  across  a  considerable  energy  range.  These  are 
the  two  sulfur  lone  pairs;  through  them  electrons  are  donated  to 
the  metal  surface.  On  the  other  hand,  the  weaker  interaction  of 
the  next  two  a"  levels  (3)  leads  to  weakening  of  the  S-C  bond, 
especially  through  the  lower  level,  which  has  strong  S-C  bonding 
character.  Although  the  band  derived  from  the  LUMO  gets  much 
wider,  there  is  only  a  very  small  electronic  transfer  from  the 
metal  to  the  adsorbate.  The  net  flow  of  electrons  occurs  in  the 
opposite  direction. 


(dotted  line)  of  a  layer  of  adsorbate  before 
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Reactions  on  the  surface 

We  saw  in  the  preceding  sections  the  extent  to  which 
ethylene  sulfide  and  trimethylene  sulfide  were  activated  by 
interaction  with  the  Mo (110)  surface.  In  ethylene  sulfide, 
activation  towards  formation  of  free  ethylene  and  sulfided 
molybdenum  was  suggested  by  the  overlap  population  changes, 
indicative  of  S-C  bond  weakening  and  C-C  bond  strengthening. 
However,  nothing  comparable  was  found  for  trimethylene  sulfide. 
Although  adsorption  may  induce  S-C  bond  weakening,  the  two  carbon 
atoms  that  are  to  form  the  new  C-C  bond  are  too  far  apart  and 
there  is  only  a  very  small  and  negative  overlap  population 
between  them.  This  may  indicate  that  the  reaction  of  adsorbed 
ethylene  sulfide  will  be  easier. 

We  decided,  however,  to  try  and  look  for  a  potential  energy 
surface  for  the  two  reactions,  in  order  to  be  able  to  compare 
their  characteristics.  The  same  reaction  as  the  one  we  used  with 
the  discrete  molecular  models  was  assumed.  There  are,  of  course, 
the  usual  limitations  involved  in  using  the  extended  Hiickel  type 
of  calculations  to  gain  information  on  how  a  reaction  occurs. 

Figure  9  gives  the  change  in  total  energy  along  the 
reaction,  and  the  variation  in  S-C,  C-C  and  Mo-S  overlap 
populations  for  the  reactions  of  both  ethylene  sulfide  and 
trimethylene  sulfide.  The  first  thing  to  notice  is  how  similar 

Figure  9  here 

to  16  and  Figures  2,  3  and  4  these  curves  are,  suggesting 
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analogous  evolutions  of  the  different  levels,  even  if  there  are 
so  many  more  for  a  surface.  As  we  saw  before,  in  the  Walsh 
diagram  of  Figure  2,  the  level  that  will  end  as  ethylene  n  is 
greatly  stabilized  and  drops  below  the  d-block  levels  just  after 
the  transition  state.  We  can  detect  this  happening  by  the 
increase  of  the  C-C  overlap  populations  after  this  level  becomes 
filled. 

The  same  happens  for  the  reaction  on  the  surface  (Figure 
9a) .  We  can  then  imagine  that  an  orbital  that  will  be  ethylene  n 
is  going  to  drop  below  the  Fermi  level  at  some  state  of  the 
reaction.  It  is  then  occupied,  and  the  C-C  overlap  population 
increases,  to  become  a  characteristic  one  for  a  C-C  double  bond. 
Unfortunately,  we  do  not  have  any  easy  way  of  constructing 
something  analogous  to  a  Walsh  diagram  for  the  reaction  on  the 
surface. 

We  tried,  however,  to  look  at  one  of  the  orbitals  of  the 
adsorbate  that  appears  to  play  a  crucial  role  in  the  reaction: 
the  ax  level  that  crosses  other  metal  d  ax  orbitals  and  finally 
becomes  the  bonding  n  level  of  ethylene  (Figures  2  and  3) .  We 
see  a  projection  of  that  M.O.  at  various  stages  of  the  reaction 
(given  by  the  stretching  of  the  S-C  bond,  as  previously,  in 
Figure  10;  note  that  it  is  not  a  linear  scale),  but  we  omitted 
the  initial  state.  The  chosen  M.O.  is  outside  the  energy  window, 

Figure  10  here 

due  to  its  strongly  S-C  and  C-H  antibonding  character. 


Energy  (eV) 


0.47 
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W 

Figure  10.  Projection  of  the  a:  molecular  orbital  of  ethylene 

sulfide  that  becomes  ethylene  ir  after  reaction  on  the 
Mo(110)  surface,  at  different  stages  of  the  reaction 
measured  by  the  stretching  of  the  S-C  bond.  The  total 
densities  of  states  (surface  and  adsorbate)  are  shown 
at  the  two  extremes  considered. 
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Stretching  the  S-C  bond  and  moving  the  hydrogens  towards  the 
nodal  plane  of  the  C  p  orbital  (see  12)  makes  it  more  and  more 
stable  as  the  reaction  proceeds.  Let  us  imagine  that  we  consider 
the  median  of  each  projected  peak  and  plot  that  energy  as  a 
function  of  the  S-C  bond  stretching.  The  resulting  curve  is 
similar  to  what  would  happen  in  the  Walsh  diagrams  of  Figure  2 
and  3  in  the  absence  of  levels  crossings.  In  the  beginning  of 
the  reaction,  the  energy  of  this  orbital  is  too  high  for  it  to 
interact  much  with  the  surface.  As  the  reaction  occurs,  we 
observe  a  spread  in  energy,  indicative  of  a  stronger  interaction. 
In  the  end  that  dispersion  disappears,  as  the  M.O.  becomes  n  and 
the  ethylene  molecule  is  very  far  away  from  the  surface. 

Let  us  now  compare  the  two  reactions.  The  decomposition  of 
ethylene  sulfide  occurs  almost  without  any  activation  energy.  On 
the  other  hand,  there  is  an  energy  barrier,  although  small,  for 
the  conversion  of  trimethylene  sulfide  to  cyclopropane.  This 
agrees  extremely  well  with  the  estimated  barriers,  based  on 
experimental  results  of  <6  and  13  kcal/mole,  respectively.3* 

Concerning  the  evolution  of  the  bonds,  we  see  that  as  the  S- 
C  bond  is  broken,  the  C-C  bond  is  strengthened  or  formed.  Also, 
and  this  is  slightly  different  for  the  surface,  two  of  the  Mo-S 
bonds  become  stronger  when  the  S-C  bonds  break.  This  effect  may 
be  an  important  one  in  contributing  to  lowering  the  activation 
energy  for  the  reaction. 

Experimental  results  show  that  the  evolution  of  ethylene  is 
slowed  down  as  the  reaction  proceeds  and  the  molybdenum  surface 
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is  getting  more  and  more  sulfided.  In  order  to  test  this,  we 
repeated  the  calculations  for  the  reaction,  adding  to  the  surface 
sulfur  atoms  in  a  1/4  coverage.  The  sulfur  atoms  were  made  to 
occupy  3-fold  sites,  as  that  is  their  usual  preference.  Both 
energy  and  overlap  population  change  as  they  did  in  the  absence 
of  added  sulfur.  If  these  results  are  to  be  trusted,  it  means 
that  the  slowing  down  of  the  reaction  is  due  not  to  the 
electronic  influence  of  co-adsorbed  sulfur,  but  to  the  lack  of 
suitable  sites  for  the  ethylene  sulfide  molecules  to  adsorb  and 
react . 

Let  us  now  turn  to  the  possibility  of  other,  non-concerted 
reactions  occurring  on  the  adsorbed  cyclic  sulfides.  In  these 
one  imagines  only  one  of  the  S-C  bonds  is  broken  at  first, 
instead  of  two.  Such  an  asymmetry  may  be  induced  by  hydrogen. 

The  presence  of  hydrogen  does  not  change  the  kinetics  of  the 
concerted  reaction.  However,  one  of  the  alternative  paths,  B  in 
1,  requires  hydrogen  in  order  that  a  terminal  methyl  group  be 
formed.  We  wanted  to  see  if  hydrogen  coadsorption  also 
influences  path  A.  We  note,  however,  that  experimentally30 
coadsorption  of  atomic  hydrogen  on  trimethylene  sulfide  does  not 
alter  the  kinetic  formation  or  the  yield  of  cyclopropane, 
suggesting  that  hydrogen  assisted  C-S  bond  activation  is  not 
occurring.  In  contrast  the  presence  of  surface  hydrogen 
increases  the  kinetics  and  yield  of  propane  formed  via  the 
thiolate  intermediate,  suggesting  the  possibility  of  hydrogen- 
assisted  C-S  bond  weakening  in  pathway  B.  We  introduced  the  one 
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hydrogen  atom  per  unit  cell  (coverage  1/4)  in  the  two  following 
situations:  adsorbate  in  3-fold  site,  perpendicular,  with  C-S-C 

plane  parallel  to  xz  and  hydrogen  in  another  3-fold  site,  20; 
adsorbate  in  3 -fold  site  but  bent  50*  towards  y  from  its 
position,  with  hydrogen  on  top  of  Mo2,  21. 


In  20,  the  hydrogen  coadsorbate  is  quite  removed  from  the 
trimethylene  sulfide  molecules.  There  are  no  steric  repulsions. 
But  on  the  other  hand,  there  is  no  activation  either.  The 
intramolecular  and  Mo-S  overlap  populations  do  not  change  when 
hydrogen  is  coadsorbed.  The  situation  changes,  however,  for  21. 
Now  the  hydrogen  is  very  close  to  one  of  the  CH2  groups  and  this 
arrangement  is  energetically  unfavorable.  The  asymmetry 
introduced  by  the  hydrogen  produces  one  slightly  weaker  S-C  bond, 
with  overlap  population  0.611,  and  a  stronger  one  (overlap 
population  0.637),  compared  to  the  free  adsorbate.  Also,  a 
strong  bond  between  the  hydrogen  and  the  carbon  of  the  nearest 
methylene  group  develops  (overlap  population  0.332).  This  is  not 
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surprising  given  the  short  C-H  distance,  1.189A.  We  can 
tentatively  assume  from  this  that  if  there  is  hydrogen  on  the 
surface  and  if  it  approaches  the  adsorbate  in  the  proper  manner 
the  formation  of  a  C-H  bond  induces  breaking  of  preferentially 
one  S-C  bond. 

Another  way  of  obtaining  the  same  result  is  by  considering  a 
tilting  movement  of  the  adsorbed  sulfide,  as  sketched  in  22, 
where  the  ring  on  the  right  is  tilted  by  6. 


22 


When  we  do  this,  the  energy  increases  guite  slowly  with  6, 
especially  for  ethylene  sulfide.  One  of  the  S-C  bonds  weakens 
with  the  simultaneous  strengthening  of  the  other  one.  This  was 
calculated  for  the  sulfur  ligand  in  a  3-fold  site,  the  ring  plane 
lying  parallel  to  xz  (18e) .  As  the  ring  is  tilted,  a  carbon  atom 
approaches  one  of  the  Mo  atoms,  23. 
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This  approach  may  facilitate  transfer  of  a  hydrogen  to 
carbon  with  breaking  of  the  S-C  bond  and  formation  of  a  thiol ate 
intermediate . 

From  this,  the  process  seems  more  likely  for  the  three- 
membered  ring,  where  we  know  that  it  does  not  place.  A  possible 
explanation  is  that  the  competing  reaction  to  form  ethylene  by  a 
concerted  process  occurs  more  easily  (there  was  no  barrier  for 
it) .  However,  for  trimethylene  sulfide,  there  is  some  barrier  to 
form  cyclopropane  concertedly.  A  competing  mechanism  then  has  a 
better  chance. 
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Discussion 

Our  calculations  suggest  that  cyclic  sulfides  prefer  to 
adsorb  on  2-  or  3-fold  sites,  a  distinction  between  these  being 
impossible  to  make  on  the  basis  of  our  results.  Similar  bonding 
preferences  are  expected  for  longer  alkyl  chain  sulfides,  as  the 
nature  of  the  frontier  orbitals  of  the  adsorbates  will  not  be 
greatly  changed.  The  bonding  takes  place  essentially  by  donation 
of  electrons  from  sulfur  lone  pairs  to  molybdenum  empty  levels. 
Only  in  ethylene  sulfide  is  there  a  significant  contribution  to 
the  bonding  from  a  S-C  bonding,  C-C  antibonding  orbital,  which  we 
think  is  responsible  for  activating  the  molecule  towards  ethylene 
formation.  We  found  that  back  donation  from  the  Mo (110)  surface 
to  the  LUMO  of  the  cyclic  sulfides  plays  only  a  minor  role  in  the 
bonding. 

The  interaction  of  cyclic  sulfides  with  a  different  surface, 
Cu (110) 17 ,  was  studied  by  means  of  ultraviolet  photoelectron 
spectroscopy,  and  it  differs  markedly  from  that  of  Mo (110) ,  in 
that  only  ethylene  sulfide  reacts  to  give  off  ethylene  (which  was 
not  detected)  and  surface  sulfur,  while  the  other  sulfides  just 
adsorb.  The  authors  describe  the  bonding  qualitatively  as  being 
mainly  electron  donation  from  the  sulfur  lone  pair  to  empty  4s, 

4p  orbitals  on  copper,  accompanied  by  back-donation  from  copper 
to  adsorbate  states  that  fell  below  the  Fermi  level.  They  do  not 
try  to  make  a  correlation  with  the  reactivity. 

The  ethylene  sulfide  ring,  eventually  containing 
substituents,  is  very  easily  desulfurized  by  reagents  such  as  RLi 
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and  Fe2(CO)g,  Fe3(CO)1218  or  RP=W(CO)5,19  for  instance,  where  the 
ethylene  that  is  formed  reacts  with  the  phosphorus  atom.  This 
reactivity  agrees  with  the  results  we  obtained  for  the  molybdenum 
molecular  models,  although  we  must  keep  in  mind  the  existing 
differences. 

A  reaction  that  comes  closer  to  the  one  on  the  surface 
involves  cyclic  sulfides  and  molybdenum  vapor.20  Ethylene 
sulfide  is  cleaved,  with  evolution  of  gaseous  ethylene  and 
formation  of  a  molybdenum  sulfide  species.  However,  as  opposed 
to  the  reaction  on  Mo (110) ,  trimethylene  sulfide  gives  only 
cyclopropane  and  propene.  The  mechanism  of  the  reaction  is  far 
from  well  understood,  but  apparently  a  biradical  intermediate  is 
not  involved. 

Another  interesting  aspect  is  the  suggestion  of  the 
formation  of  a  pentacoordinate  sulfur  intermediate  preceding 
olefin  formation  when  ethylene  sulfides  react  with  RLi . 18  This 
is  comparable  to  the  coordination  of  sulfur  on  the  3-fold  site 
that  we  tentatively  propose  to  be  one  of  the  most  activating. 

Much  work  has  been  done  in  order  to  understand  how  reactions 
occur  on  surfaces  and  why  several  crystal  faces  exhibit  different 
adsorption  and  reactivity  patterns,  an  example  being  the 
reduction  of  nitric  oxide,  NO,  on  platinum.21  A  few  models  have 
been  proposed,  which  deal  with  the  characteristics  of  the 
adsorbate-surface  interaction  in  the  initial  stages  of  the 
reaction.  On  the  other  hand,  it  was  suggested,  although  in  the 
context  of  another  type  of  catalytic  reactions,  that  the 
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interaction  with  the  transition  state,  not  with  the  initial 
state,  is  determinant  of  the  catalyst  activity.22  This  idea  was 
not  very  straightforward  to  apply  to  the  comparison  of  the 
reactions  of  ethylene  and  cyclopropane  formation  from  CZH*S  and 
C3H6S,  respectively,  on  Mo (110) .  One  of  the  difficulties,  and 
this  is  different  from  other  studied  reactions, 15b,n  is  that  the 
gaseous  ethylene  products  formed  are  not  trapped  by  the 
surface. 3b,e,f  The  surface  only  acts  on  the  reaction  through  the 
sulfur  atom  of  the  cyclic  sulfide.  This  was  the  reason  that  led 
us  to  follow  a  specific  pathway  for  the  two  reactions,  even  with 
all  the  limitations  involved.  As  the  ring  size  increases,  more 
severe  geometric  changes  are  needed  for  a  concerted  reaction, 
increasing  the  activation  energy,  thus  making  that  reaction  path 
less  and  less  competitive  with,  for  instance,  selective  breaking 
of  only  one  S-C  bond.  We  saw  earlier  how  this  could  be  induced 
either  by  coadsorbed  hydrogen  or  by  tilting  movements  of  the 
adsorbate  which  do  not  require  too  much  energy.  A  more  detailed 
study  of  this  alternative  mechanism  is,  unfortunately,  outside 
our  calculational  possibilities,  due  to  the  enormous  number  of 
degrees  of  freedom  involved. 
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Appendix 

All  the  calculations  were  of  the  extended  Huckel  type,7  with 
the  tight-binding  approach8  for  the  surface  calculations.  The 
parameters  used  are  collected  in  Table  4. 

Table  4  here 

The  distances  taken  in  the  molecular  models  MoH*  (SCnHj  *'  and 
Mo2H6(SC2HJ 6‘  were,  respectively,  Mo-H  1.80,  Mo-S  2.44,  Mo-Mo 
2.12k.  Ethylene  sulfide  and  trimethylene  sulfide  were  considered 
to  have  the  same  geometry  as  determined  for  the  free  molecules. 
The  two-dimensional  molybdenum  surface  consisted  of  three  layers. 
The  adsorbate  coverage  was  kept  at  1/4  throughout  all  the 
calculations. 

Sets  of  6  or  8  k-points  in  the  irreducible  wedge  of  the 
oblique  unit  cell  in  reciprocal  space  were  used,  according  to  the 
symmetry  of  the  problem.  They  were  chosen  following  the 
geometrical  method  of  Ramirez  and  Bohm.23 


Table  4 


Parameters 

used  in  the 

extended 

Hiickel 

calculations 

Atom 

Orbital 

Hii(eV) 

f  2 

Cia 

c2a 

Mo 

5s 

-8.34 

1.96 

5p 

-5.24 

1.90 

4d 

-10.50 

4 . 54 

1.90 

0.5899 

0 . 5899 

s 

3s 

-20.0 

1.817 

3P 

-13 . 3 

1.817 

C 

2s 

-21.4 

1.625 

2p 

-11.4 

1.625 

H 

is 

-13.6 

1.300 

aCoef ficients  used  in  the  double  f  expansion  of  the  4d  orbitals 
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